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By comparing photoemission spectroscopy with a non-perturbative dynamical mean field theory 
extension to many-body ab initio calculations, we show in the prominent case of pentacene crystals 
that an excellent agreement with experiment for the bandwidth, dispersion and lifetime of the hole 
carrier bands can be achieved in organic semiconductors provided that one properly accounts for the 
coupling to molecular vibrational modes and the presence of disorder. Our findings rationalize the 
growing experimental evidence that even the best band structure theories based on a many-body 
treatment of electronic interactions cannot reproduce the experimental photoemission data in this 
important class of materials. 



Introduction. — Organic semiconductors are key ma- 
terials for future applications ranging from flexible elec- 
tronics to photovoltaics fl] . Because of the weakness of 
van der Waals inter-molecular bonds, these materials are 
commonly believed to stand in an intricate region be- 
tween the molecular limit and the extended band pic- 
ture, calling for concepts that go beyond the standard 
paradigms of inorganic semiconductors. Experimental 
charge transport and optical studies over the past few 
years have indeed given contrasting indications on the 
nature of the electronic carriers. The existence of low- 
mass quasiparticles indirectly inferred from optical sum 
rules ^ , and the apparently "band- like" temperature de- 
pendent electronic mobilities observed in the best crys- 
talline systems [3H5] are difficult to reconcile — within 
the framework of band theory alone — with absolute 
values of the mobility that hardly exceed few tens of 
cm? /Vs. The coupling to molecular vibrations as well 
as the presence of disorder are often regarded as the pri- 
mary causes of the poor conductive properties of this 
broad class of materials 7 , as both phenomena can slow 
down the electron motion in the already narrow bands 
constructed from vr-intermolecular overlaps. Still, no di- 
rect proof on how such microscopic mechanisms affect 
the electronic states has been given to date. 

As of today, while the strongest efforts are devoted 
to understanding and improving their charge transport 
characteristics, a proper description of even the most ba- 
sic electronic properties such as the band dispersion and 
carrier lifetime — the foundations of our understanding 
of conventional semiconductors — remains a challenge in 
organic semiconductors. Recent angle resolved photoe- 
mission (ARPES) experiments performed on crystalline 
samples [8Ul3] have evidenced that the best band struc- 
ture theories based on a many-body treatment of elec- 



tronic interactions cannot reproduce the experimental 
spectra, as they predict electronic bandwidths that are 
significantly smaller than measured experimentally. Here 
we provide a combined theoretical/experimental analy- 
sis of the ARPES spectra of pentacene — a material 
that is often regarded as a model compound for organic 
molecular solids. By focusing on the full momentum de- 
pendence of the spectral features, we demonstrate that 
the interaction of electrons with the internal vibrations 
of the molecules deeply modifies the nature of the ex- 
tended electronic states. The multiple vibrational over- 
tones that are commonly observed in angle-integrated 
photoemission spectra |14H16j translate into a complex 
fine-structure in momentum space, unveiling an unex- 
pectedly large separation of the electronic bands that is 
ultimately responsible for the observed increase of the 
total bandwidth. Analogous molecular fingerprints are 
predicted to affect the electronic bands of broad classes 
of molecular systems including organic semiconductors, 
charge-transfer organic salts and doped organic conduc- 
tors [TT]. 

Photoemission experiment. — In order to probe the ef- 
fect of electron-molecular vibration interactions on the 
electronic spectra of organic solids, ARPES experiments 
were performed on in situ grown pentacene films of thick- 
ness d > lOOA deposited on Bi(OOl). This is in contrast 
with previous experiments performed on organic mono- 
layers [HllTllIlB], where the energy- momentum dispersion 
could be seriously affected by the binding to the sub- 
strate. All experiments were carried out at the University 
of Wisconsin Synchrotron Radiation Center (SRC). A de- 
tailed description of the sample preparation process and 
of the resulting crystalline film including lattice constants 
a and b are reported elsewhere [ini HI] . The combined 
photon and electron energy resolution was Ai? ~ 40 
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FIG. 1: (a) Color density plot of experimental ARPES data. The second derivative of the data has been taken to best 
visualize the energy positions of the dispersive features. The green dots are the peak positions as derived from gaussian fits 
of the individual spectra. The blue dot corresponds to a scan with a different incoming photon energy, which recovers the H2 
spectral weight around F that is hidden due to the matrix element dependence on the photon incoming energy 18.. The gray 
line is the ab initio band dispersion from Ref. 1191 (b) Calculated spectrum in the presence of EMV interactions and disorder 
(see text). The dots are the experimental peak positions from panel (a); (c) Theoretical HOMO band dispersion, as obtained 
from the poles of the spectral function. Ab initio (gray), disordered (blue: A = 75 meV, Ep = 0) and with EMV interactions 
(red: A = 75 meV, Ep = 69 meV). The inset shows the corresponding path in the Brillouin zone. 



meV and, for 15 eV photons (used in most experiments) 
the electron momentum resolution was Afc|| < 0.03 A-i. 
Fig. [T^ shows the photoemission spectra measured at 
T — 75K near the top of the highest occupied molecular 
orbital (HOMO) as a function of parallel momentum fc|| 
along the path illustrated in the inset of Fig. [TJ;. Since 
the lattice structure has two molecules per unit cell, two 
dispersive branches can be identified, that we denote as 
HI and H2. The dots in Fig. [l^ are the positions of 
the two most prominent peaks at each given fcy , as ob- 
tained from gaussian fits to the individual spectra. The 
solid lines are the band dispersions obtained from state- 
of-the-art ab initio calculations for the crystalline film 
phase which corresponds to the pentacene polymorph of 
our sample [19, 20 . It is apparent that the separation 
between the HI and H2 branches in the experiment is 
much larger than predicted. As a consequence the mea- 
sured total HOMO bandwidth, Wexp = 450 ± 15 meV, 
defined as the distance between the HI and H2 peaks 
at the M point [T^, is also much larger than the calcu- 
lated values for the pentacene crystalline film phase, both 
within DFT ^ (Wq = 348 meV) and within the more 
accurate GW method [5T1 (Wq = 360 meV). The discrep- 
ancy is well beyond the experimental uncertainty, which 
confirms the systematic findings reported in the most re- 
cent ARPES measurements on pentacene and rubrene 
[8HT5]. This indicates that the band-structure calcula- 
tions alone, even including electronic many-body corre- 
lations as is done in Ref. [2TJ are not sufficient to properly 
describe the observed excitation spectra. 



Theoretical modeling. — In order to solve the discrep- 
ancy between theory and experiment, we now calcu- 
late the photoemission spectral function by consider- 
ing the coupling of holes in the HOMO band with the 
high-frequency, intramolecular vibrations arising from 
the carbon-carbon stretching forces. We use a non- 
perturbative, theoretical method in order to avoid any a 
priori assumption on the relative importance of the dif- 
ferent phenomena at work. This is a crucial issue since in 
organic semiconductors the characteristic energy scales, 
namely (a) the electronic bandwith Wq, proportional to 
the intermolecular electronic transfer energies, (b) the 
strength of the electron-EMV coupling, measured by the 
structural relaxation energy Ep associated with adding a 
charge to a molecule, (c) the molecular vibrational ener- 
gies and (d) the energy fiuctuations due to disorder, that 
cannot be ignored when dealing with real systems, are all 
of comparable magnitude (^ 0.1 eV). This means that no 
small parameter can be identified that would allow appli- 
cation of perturbation theory. The microscopic parame- 
ters of the EMV interaction specific to pentacene are ob- 
tained here via accurate many-body ab initio calculations 
within the GW formafism [fgll^^tfcs] . This method has 
been shown to yield much better electron-phonon cou- 
pling potentials as compared to standard DFT-LDA cal- 
culations, in the case of 7r-conjugated fullerenes [24! and 
graphene [53]. We then calculate the resulting photoe- 
mission spectrum via DMFT [211 HZ] , which is able to ad- 
dress in an unbiased way the intricate physical regime of 
interest in organic semiconductors. The DMFT scheme 
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also allows one to include the presence of disorder that 
is expected to be relevant in experimental samples, in 
the form of spatial fluctuations of the molecular energy 
levels. Experimental estimates of the energy spread A 
range from few tens of meV in single crystals to A ~ 0.1 
eV or above in amorphous samples [55]. Our crystalline 
films should be located between these two limits as in- 
dicated by the sharpness of our reflection high energy 
electron diffraction (RHEED) patterns. Finally, the in- 
teraction with low-frequency vibrations of the molecules 
[Tll[53|, in the range hw < 40 meV, is effectively included 
in the calculation via an increased value of A. These low- 
frequency modes related both to the large molecular size 
and to the mechanical softness of the material do not give 
rise to observable features in the photoemission spectra, 
as their energy lies below the experimental resolution, 
but they do act as an additional intrinsic source of fluc- 
tuation for the molecular energy levels [30l |3T] . 

Molecular origin of the band separation. — Fig. [T]d is a 
color density plot of the hole spectral function obtained 
from the theory including both molecular and itinerant 
aspects on the same footing, as described in the preced- 
ing paragraph. We take the calculated value Ep = 69 
meV for the molecular relaxation energy (corresponding 
to a reorganization energy 2Ep = 138 meV), hil = 174 
meV for the average energy of intra-molecular vibrations 
[TH], a non-interacting bandwidth Wq = 348 meV and 
assume a total disorder strength A = 75 meV. The sim- 
ilarity with the experimental spectra is remarkable. The 
most striking point is that the separation between the HI 
and H2 bands is much larger than in the band-structure 
prediction, restoring the agreement with the locus of the 
experimental peaks (dots). Also in agreement with the 
experiment, the features of the calculated spectrum ex- 
hibit a considerable broadening, a point that will be dis- 
cussed below. 

To ascertain the microscopic origin of the observed 
H1/H2 separation. Fig. ^ shows the HOMO band dis- 
persion obtained from the theory upon the subsequent 
inclusion of disorder (blue curve) and EMV interactions 
(red curve) on top of the band structure as parametrized 
in Ref. [1^1 Disorder in the molecular energy levels causes 
a rather uniform increase of the band dispersion, but does 
not modify the overall band shape [32]. Instead, adding 
the interaction with high-frequency molecular vibrations 
causes a sizable upwards shift of the HI band, leaving 
the H2 band dispersion essentially unchanged from the 
non-interacting case. This shift is caused by the molec- 
ular relaxation energy gain Ep, that is reflected in the 
extended states of the solid through a stabilization of the 
HI branch. We conclude that the interplay between EMV 
coupling and disorder is at the origin of the large H1/H2 
separation observed in the experiment. The present anal- 
ysis confirms the results obtained in Ref. [3^ based on a 
related one-dimensional model. 

Spectral hallmarks of EMV interactions. — Having 
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FIG. 2: (a) Photoemission spectrum measured at the M 
point. The dashed line is a fit of the data with three gaus- 
sian peaks, (b) Calculated spectrum and fit. Parameters are 
Ep — 69 meV and A = 75 meV. W is the total bandwidth 
defined as the distance between the HI and H2 peak, (c) The 
fine structure induced by the EMV interaction is unveiled by 
removing the disorder in the calculation, (d) The calculated 
spectrum in the absence of EMV interactions, but with the 
same degree of disorder. In (b,c,d) the spectrum has been 
convoluted with a Gaussian of FWHM = 40 meV to mimic 
the experimental resolution. 



shown that the main dispersive features of the ARPES 
data can be explained by the proper inclusion of inter- 
actions with high-frequency molecular vibrations and a 
moderate amount of disorder, we now proceed to show 
that additional hallmarks of the EMV coupling are seen 
in the angle-resolved spectra, providing support to the 
proposed scenario. Fig. [2^ reports the photoemission 
intensity measured at the M point (full line), and Fig. 
[2]d is the corresponding calculated spectrum. In addition 
to the HI and H2 peaks that are clearly resolved, an ex- 
tra feature can be actually recognized, whose existence is 
substantiated by fitting the spectra with three gaussian 
peaks (dashed lines). Being located at a distance ~ Ti^ 
from the main HI peak, it is tempting to associate this 
feature with a vibrational overtone of the main electronic 
excitation. That this is indeed the case is demonstrated 
by repeating the theoretical calculation in the absence of 
disorder, which reveals all the fine structure that is oth- 
erwise smeared out by the energy fluctuations (Fig. [2j;). 
The overtone actually disappears as expected when the 
EMV interaction is turned off (Fig. [2ji). 
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FIG. 3: (a) Calculated HOMO bandwidth as a function of 
the EMV coupling strength for different values of the disorder 
parameter. The shaded area is the measured value (see text). 
The horizontal bar indicates the range of calculated Ep values 
in the literature, (b) Renormalization of the HI band along 
Pa2 (M — f distance, see Fig. [T]). The dashed line is the 
polaronic renormalization expected in the molecular limit. 



Our analysis shows that the presence of multiple vibra- 
tional overtones is also at the origin of the large broad- 
ening of the H2 peaks observed in the experiment, whose 
tails extend to very large binding energies (see both Fig. 
[2^ and [2]d) . The reduced lifetime of the H2 states again 
results from a non-trivial interplay between vibrational 
shakeoff excitations and disorder. The individual reso- 
nances contributing to the H2 weight are not resolved in 
the experiment, as these merge into a single broad feature 
due to the presence of disorder. Such vibrational over- 
tones are however clearly revealed in the calculation on 
a perfectly ordered crystal (Fig. [2]:) and could possibly 
be detected in future experiments on cleaner samples. A 
more detailed discussion on the interplay between EMV 
interactions and disorder, not limited to the representa- 
tive M point considered here, is presented in [18] . 

Quantitative analysis and discussion. — Finally we 
show that combining model Hamiltonian studies with ab 
initio calculations as we have done here can actually be 
used to extract quantitative information on the relevant 
microscopic parameters of the materials directly from the 
experiment. Fig. [3^ illustrates the evolution of the to- 
tal calculated HOMO bandwidth W resulting from the 
molecular relaxation phenomenon as a function of the 
EMV coupling strength A = 2Ep/Wo, for different val- 
ues of the disorder parameter A. Taking the calculated 
value, Ep ~ 69 nieV, our experimental data are compat- 
ible with a total disorder strength A = 75 ± 15 meV, in 
agreement with the values available in the literature [55] . 
As mentioned previously, this is an effective parameter 
that also includes the intrinsic fluctuations originating 
from the low- frequency molecular vibrations, whose con- 
tribution can be estimated to ~ 20 meV [211 [50] ■ The 
estimate A = 75 meV given above therefore constitutes 
an upper bound to the actual structural disorder present 



in the sample. We note that the value of the molecular 
relaxation energy itself could be extracted from the ex- 
periment if a precise independent estimate of the disorder 
fluctuations were available. 

Beyond these quantitative considerations, the present 
study clearly establishes the theoretical framework that 
should be worked with in upcoming studies of transport 
properties, since a good description of the momentum 
and energy dependence of carriers is a necessary basis 
to properly assess mobility properties. In particular, our 
results show unequivocally that the total HOMO band- 
width increases with the EMV coupling. This behavior 
is opposite to the bandwidth decrease that constitutes 
the basis of most theories of charge transport in organic 
semiconductors, but that is only expected to apply in the 
molecular limit, Wq <C HQ [3 [331 [SI]. Upon closer inspec- 
tion, we find that a partial band shrinking occurs in a re- 
stricted energy range, roughly within a window hfl from 
the top of the band (the ground state for holes) . The 
width of the HI branch along Pa2, reported in Fig. [Sjj, 
indeed shows a moderate decrease upon increasing the 
EMV coupling, which however does not obey the text- 
book |34j exponential renormalization cx exjp{—Ep/Hn) 
(dashed line). The fact that the phenomenon of band 
narrowing is shown here to play a reduced role should be 
taken into account when constructing a theory for charge 
transport in organic semiconductors. 

Concluding remarks. — There is growing experimen- 
tal evidence that the electronic properties of organic 
molecular semiconductors cannot be understood within 
the commonly accepted band structure picture that pre- 
vails for inorganic semiconductors. This discrepancy en- 
tails that a consistent description of the extended states 
in molecular solids should take full account of the mi- 
croscopic processes taking place at the molecular level. 
The analysis presented here indeed demonstrates that 
the proper inclusion of the interaction with molecular 
vibrations and disorder, beyond electronic band theory 
calculations, provides a remarkably accurate description 
of the experimental photoemission spectra in pentacene. 
The coexistence of dispersive features characteristic of 
the extended band regime, together with multi-phonon 
shakeoff resonances reminiscent of the molecular spectra, 
provides solid spectroscopic evidence for the widespread 
idea that organic materials are located in a crossover re- 
gion where the band and molecular characters are inex- 
tricably linked. 
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